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A quasi-particle description of (2 + 1)- flavor lattice QCD equation of state 
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A quasi-particle model has been employed to describe the (2 + l)-flavor lattice QCD equation of 
state with physical quark masses. The interaction part of the equation of state has been mapped 
to the effective fugacities of otherwise non-interacting quasi-gluons and quasi-quarks. The mapping 
is found to be exact for the equation of state. The model leads to non-trivial dispersion relations 
for quasi-partons. The dispersion relations, effective quasi-particle number densities, and trace 
anomaly have been investigated employing the model. A Virial expansion for the EOS has further 
been obtained to investigate the role of interactions in quark-gluon plasma (QGP). Finally, Debye 
screening in QGP has been studied employing the model. 
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I. INTRODUCTION 

The main purpose of this article is to explore the quasi- 
particle picture of Quantum-Chromodynamics (QCD) at 
high temperature. In particular, we wish to describe re- 
cent lattice data on (2+l)-flavor QCD equation of state 
(EOS)by employing a quasi-particle model The 
EOS is an important quantity to study the properties of 
hot QCD matter which is commonly known as quark- 
gluon plasma (QGP) in relativistic heavy ion collisions 
at BNL, RHIC and CERN, LHC. This study is man- 
dated by the strongly interacting nature of QGP which 
has been inferred from the recent experimental observa- 
tions at RHIC 0-0| • This observation is consistent with 
the lattice simulations of the EOS which predict a 

strongly interacting behavior even at temperatures which 
are a few T c (QCD transition temperature). 

The most striking features of RHIC results Q are the 
large collective flow and strong jet quenching of high 
transverse momentum jets shown by QGP. Similar con- 
clusions have been drawn from the very recent prelimi- 
nary results for Pb-Pb collisions at LHC at s/s = 2.76 
TeV fl5l - fl8l ]. This has led to a tiny value of the shear 
viscosity to entropy density ratio for QGP and near per- 
fect fluid picture of QGP [a, EH2I (except near the QCD 
transition temperature where the bulk viscosity of QGP 
is equally important as shear viscosity (2714291 ] ) . In an at- 
tempt to appreciate this result, interesting an alog ies have 
been drawn with ADS/CFT correspondence [30( as well 
as with some strongly coupled classical systems [3l| . In 
any case, the emergence of strongly interacting behavior 
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puts into doubt the reliability of a large body of analyses 
which are based on ideal or nearly ideal behavior of QGP 
in heavy ion collisions. 

In the light of these observations, it would be right 
to state that QGP may lie in the strongly interacting 
domain (non-perturbative) of QCD. Therefore, lattice 
gauge theory [32H34| would be the best approach to ad- 
dress the physics of QGP in RHIC in terms of a reliable 
EOS which is very precisely evaluated The lattice 

EOS is far from being ideal. The EOS is rj 10% away 
from its ideal behavior even at 4T C . However, in sev- 
eral works devoted to QGP ideal EOS is employed. This 
is certainly not desirable for QGP in RHIC. Therefore, 
there is an urgent need to address this issue by develop- 
ing models to employ realistic QGP equations of state 
to investigate the bulk and transport properties of QGP. 
This could be done by casting hot QCD medium effects 
in terms of effective quasi-particle degrees of freedom. 

There have been several attempts to describe QCD 
medium effects at high temperature in terms of quasi- 
particle degrees of freedom. These attempts include, 
(i) the effective mass approaches to study QCD ther- 
modynamics I35H39B . and (ii) approaches based on the 
Polyakov loop pfj - l44| . A different approach in terms of 
quasi-particles, inspired by Landau theory of Fermi liq- 
uids has been proposed recently both for EOS based on 
pqcd iBiEi and pure lattice guage theory 0). 
This model is fundamentally different from the above 
two approaches, and quite powerful: it reproduces the 
EOS with remarkable accuracy, especially in the case of 
lattice EOS-where it is exact; the collective nature of 
the quasi-gluons is manifest in the single particle disper- 
sion relations. It is also successful in terms of predictions 
regarding the bulk and transport properties of QGP [l|- 
lUlUl. Refs. 0, H3 showed that the shear viscosity and 
its ratio with the entropy density (77, tj/s) arc highly sen- 
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sitive to the interactions. They could be taken of as good 
diagnostics to distinguish various EOS at RHIC. 

The model was tested only against the pure SU(3) 
gauge theory EOS [||, where the description was, as we 
mentioned, exact. It was not employed for extracting a 
quasi-particle description in the case of full QCD, by the 
inclusion of quark sector. In this paper, we remedy this 
draw back and extend the model to the matter sector, 
by taking up a recently computed (2 + l)-flavor lattice 
QCD equation of state with physical quark masses [Ioj |. 
It is to be noted that this EOS has further been refined 
by improving the accuracy in [ill . Il2j . Here, we adopt 
the philosophy same as in 0. Again, we map (2 + 1)- 
flavor lattice QCD data for the EOS [l(| in terms of 
quasi-particle degrees of freedom which are free up to 
effective gluon-fugacity, z g and effective quark fugacity, 
z q . In this model, the strange quark sector is different 
from light quark sector due to contributions coming from 
the strange-quark mass. They are otherwise characterize 
by the same effective fugacity z q . Such a characterization 
could only be possible because the mass corrections from 
light quark-sector in the deconfined phase of QCD are 
very very small. So, the model will be more realistic at 
higher temperatures. 

Here, it is worth mentioning that such a quasi-particle 
description of the lattice QCD EOS could be thought of 
as a first step towards an effective field theory/ effective 
kinetic theory to explore complicated nature of strong in- 
teraction in QGP. Leaving these ambitious investigations 
for future studies, here, we have attempted to understand 
the role of QCD interactions in terms of a Virial expan- 
sion for QGP employing the quasi-particle model. The 
Virial expansion has been obtained in terms of effective 
quasi-particle number densities which are ,in turn, ex- 
pressed in terms of z g / q . As we shall see that the Virial 
expansion of the EOS is very helpful to understand the 
role of strong interaction in QGP and may perhaps play 
crucial role in developing the effective models. 

The paper is organized as follows. In section II, we 
introduce the quasi-particle model and study its features 
and physical significance. Here, we discuss the viability 
of the model by studying the temperature dependence of 
the quasi-particle pressure, and trace anomaly in terms of 
effective fugacities. We find that the model reproduces 
the EOS almost exactly. In section III, we discuss the 
physical significance and viability of the quasi-particle 
model. In section IV, we discuss the implications of the 
model. Here, we propose a Virial expansion for QCD 
at high temperature, in terms of effective quasi-particle 
number densities and explore the role of interaction in hot 
QCD. We further study the Debye screening and charge 
renormalization in hot QCD. In section V, we present the 
conclusions and the future direction of the work. 



II. THE QUASI-PARTICLE MODEL 

Before, we introduce the quasi-particle description, let 
us define the notations. The quantities, z g , and z q will 
denote effective gluon and quark/anti-quark fugacities 
respectively. The quasi-gluon equilibrium distribution 
function will be denoted by quasi-quark/anti-quark 
distribution function by f% for light quarks (u,d), and 
fl q for strange quark. The respective dispersions (single 
quasi-particle energy) will be denoted as and w*. 

n g,q,s denotes the effective quasi-particle number densi- 
ties. In all the physical quantities that will be discussed 
in the paper, the subscript g will denote the gluonic con- 
tribution while q and s denote the contributions from the 
light quark sector and strange quark sector respectively. 

A. Underlying distribution functions and effective 
fugacity 

We initiate the model with the ansatz that the Lattice 
QCD EOS can be interpreted in terms of non-interacting 
quasi-partons having effective fugacities which encodes 
all the interaction effects. In the present case, we have 
three sector, viz., the effective gluonic sector, the light 
quark sector, and the strange quark sector. Here, the ef- 
fective gluon sector refers to the contribution of gluonic 
action to the pressure which also involves contributions 
from internal fermion lines. Due to purely phenomeno- 
logical reason, this sector can be rccastcd in terms of ef- 
fective gluon quasi-particles (which are free gluons with 
effective fugacity). Similary the other two sectors also 
involve interactions among quark, anti-quarks, as well 
as their interactions with gluons. The effective gluon 
fugacity, z g is introduced to capture the interaction in 
the effective gluonic sector. On the other hand, z q cap- 
tures interactions in other two sectors. The ansatz can 
be translated to the form of the equilibrium distribution 
functions, /f , /| , and /| g as follows, 

z g exp(-(3p) 

Jeq / \ ) 

( 1 - z g exp(-fip) J 
, g z g exp(-/3p) 

Jeq / \ J 

{l + z q exp(-#p) J 

z q exp(~(3^p 2 +TO 2 ) 

Jeq / \ J \L) 

il + z q cxp(-/V;P 2 + m 2 ) J 

where m denotes the mass of the strange quark, which 
we choose to be O.lGeV. /3 = T^ 1 denotes inverse of 
the temperature. Note that we are working in the units 
where Boltzmann constant, Kb = 1, c = 1, and h/2w = 
1. We use the notation v g = 2(iVj — 1) for gluonic degrees 
of freedom .^ = 2x2x^x2 for light quarks, v s — 
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2 x 2 x N c x 1 for the strange quark for SU(N C ). Here, 
we are dealing with ££7(3), so N c = 3. Since the model 
is valid in the deconfincd phase of QCD (beyond T c , T c 
is the QCD transition temperature), the masses of the 
light quarks can be neglected. Therefore, in our model 
we only consider the mass for the strange quarks. 

As it is well known that QCD thermodynamics at high 
temperature is described in terms of a Grand canonical 
ensemble. Now, it is straight forward to write down an 
effective Grand canonical partition function for hot QCD 
which yields the forms of the distribution function given 
in Eq.((T]). We denote the partition function by Z = (Z g x 
Z q x Z s ). The corresponding expressions in terms of z g 
and z q are as follows, 



d 3 p 
8^ 

d 3 p 
8^ 



HZ g ) = -v g V / ln(l - z g exp(-/?p)) (2) 



ln(Z q ) = Vq V \ ln(l + z q cxp(-/3p)) (3) 



ln(Z s ) = v s V I ln(l + z 9 exp(-/?(v/p 2 + m2))(4) 



ln(Z) = ln{Z g ) + ln(Zq) + ]n(Z a 



(5) 



Now using the well known thermodynamic relation, 
P/3V = hi(Z), we can match the rhs of Eq.(5), with the 
lattice data for the pressure for (2 + 1)- flavor QCD \v$, 
where P denotes the pressure and V denotes the vol- 
ume. From this relation, we can in principle determine 
the temperature dependence of z g and z q . As emphasized 
earlier, z g is determined from the contribution to the lat- 
tice pressure purely from gluonic action. This particular 
contribution to the pressure is denoted as P g . Remain- 
ing part of the pressure is utilized to fix the temperature 
dependence of z q . Now, we have two relations and two 
unknowns. Next, we discuss the determination of z g and 

Zn- 




FIG. 1. (Color online) Behavior of z g as a function of T/T c . 
The points denote the values obtained from lattice data and 
solid line denote the fitting function. The fitting parameters 
are listed in Table. I. 




1. Determination of z g and z q 

We determine z g and z q numerically. z g has been de- 
termined using the relation, 

P g = ^ 2 9 J du u 2 ln(l - z g exp(-u)). (6) 

On the other hand, z q has been determined numerically 
using the following relation, 

o-4 roo f 

( P - P a) = ^2j duu 2 lvq\n(l + z q cxp(-u)) 
+v s ln(I + z q exp(— y/ u 2 + m 2 )) |, 

(7) 

where m = /3m and u (u = f3p) is a dimensionless quan- 
tity. We have recorded those values of z g and z q which 
satisfy Eqs. ([6]) and ([7]). Next, we discuss their behavior 
with temperature. 



FIG. 2. (Color online) Behavior of z q as a function of T/T c . 
The points denote the values obtained from lattice data and 
solid line denote the fitting function. The fitting parameters 
are listed in Table. I. 



2. Behavior of z g and z q 

The determination of the the quasi-parton distribution 
functions given in Eq.([T]) is complete once the tempera- 
ture dependence of z g and z q is fixed. The behavior of z g 
and z q as a function of temperature is shown in Fig. 1 
and in Fig. 2 respectively. Clearly, both of them acquire 
their ideal values (unity) only asymptotically. At lower 
temperatures, the magnitude of both z g and z q is smaller 
indicating the larger strength of interactions there. 

For further analysis, we seek analytic forms for z g and 
z q as a function of temperature, which would render the 
computation more amenable. At this juncture, we note 
that there are limited number of lattice data for the pres- 
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TABLE I. Fitting parameters for z g and z q 



Z S,Q 


a a,i 


b 9,q 


/ 

a 9,Q 




Gluon 
Quark 


0.803±0.009 
0.810± 0.010 


1.837±0.039 
1.721±0.040 


0.978± 0.007 
0.960 ± 0.007 


0.942±0.035 
0.846 ± 0.033 



sure for a huge range of temperature ((0.5— 4.0)T C ). Since 
the effective fugacities have been determined from the 
lattice pressure, the limitations get inherited by them as 
well, which further passed on to other thermodynamic 
quantities such as energy density, entropy density and 
trace anomaly. We hope that the functional form is not 
drastically altered by future refinements in lattice data. 

We find that there is no universal functional form that 
describes the data in the full range of temperatures, ei- 
ther in the gluonic sector or the quark sector. There 
are some common features though. Both the sectors are 
characterized by a 'low temperature' and a 'high temper- 
ature' regime, with the cross over temperatures given by 
x 3t q = Tg, q /T c ~ 1.68, 1.70 respectively. The functional 
forms on cither side of the same for both the sec- 

tors, but with different parameters. Thus, when x < x qjg , 
the good fitting function has the form a gq exp(—b gjq /x 5 ). 
In the complimentary case x > x q _ g , it has the form 
a' g exp(— b' g q /x 2 ). The latter form is mandated by the 
fact that at high temperature the trace anomaly, e ~~4 P 
predominantly goes as 1/T 2 [l(j in lattice QCD. We list 
the fitting parameters for z g and z q in Table I. 

We shall utilize these forms to study temperature de- 
pendence of the trace anomaly in later part of the paper. 
We shall see that these forms correctly reproduce the 
high and low temperature behavior of the trace anomaly. 

Temperature dependences of effective fugacities, in 
Fig. 1 and Fig. 2 reveal that the effective gluon and 
quark fugacities arc of same order of magnitude for the 
whole range of temperature. This indicates that effective 
gluons and quarks contribute equally in our description. 
Its possible physical consequences, and an understand- 
ing from basic calculations in QCD is beyond the scope 
of the present work and will be a matter of future inves- 
tigations. 

It is worth noting that effective fugacity descriptions 
have been earlier employed in condensed matter systems 
in the last decade. To study the nature of Bosc-Einstein 
(BE) condensation transition in interacting Bose gases, 
a parametric EOS in terms of the effective fugacity has 
been proposed by Li et. al [if|. This provides a scheme 
to explore the quantum-statistical nature of the BEC 
transition. There have been other works to study the 
non-interacting BE systems in harmonic trap [47j as well 
interacting bosonic systems [48|. Moreover, effective fu- 
gacity description has been used for a unitary fermion 
gas by Chen et al [4!| for studying thermodynamics with 
non-Gaussian correlations. None of them employed the ef- 
fective dispersion relations which we obtain naturally in 
this work. We shall now proceed to discuss the physical 
significance of the quasi-particle model and its viability. 



III. PHYSICAL SIGNIFICANCE AND 
VIABILITY OF THE MODEL 

A. The modified dispersion relations 

It has been emphasized in Ref. that the physical 
significance of effective fugacity could be seen in terms 
of modified dispersion relations. The effective fugacities 
modify the single quasi-parton energy as follows, 

uj3=p + T 2 d T ln{z g ) 
u q p =p + T 2 d T ln{z q ) 



x/p 2 + m 2 + T 2 d T ln(z q ). 



(8) 



These dispersion relations can be interpreted as follows. 
The single quasi-parton energy not only depends upon 
the momentum but also gets contribution from the collec- 
tive excitations of the quasi-partons. The second terms 
is like the gap in the energy due to the presence of quasi- 
particle excitations. This immediately reminds us of Lan- 
dau's theory of Fermi -liquids. Therefore, it is safe to say 
that the present quasi-particle model is in the spirit of 
Landau theory of Fermi liquids. These modified disper- 
sion relations in Eq.® have emerged from the thermo- 
dynamic definition of the average energy of the system, 
due to the temperature dependent fugacities, z g / q . Let us 
consider the expression for the energy-density, e obtained 
in terms of the Grand Canonical partition function, Z as, 



1 dln(Z) 

V d/3 ' 



(9) 



Substituting for the effective partition function (Z 
we obtain, 



Zg x Z q x Z s ) 



d 3 p^p + T 2 &rH^)))f! q 
+ ^3 / d 3 p[(p + T 2 d T Hz q ))y q f t 
+ [yWT^ 2 + T 2 d T m(z 9 )) v s fl q ) 



: 3(P g + P q + P. 

T 2 T ln(z q ) 
2tt 2 



2^2 J a P "sJeq 

d 3 P Uf? q + Vsft q y (10) 



The above equation can be recasted employing the ex- 
pression for the pressure in terms of the temperature de- 
pendent, z g / q as, 



(6 - 3P) _ d ( P 

T 4 q t I T 4 



(11) 



Therefore, these modified dispersion relations natu- 
rally ensure the thermodynamic consistency condition in 
high temperature QCD, and lead to the trace anomaly, 
which we have discussed, in detail, in the next subsection. 
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FIG. 3. (Color online) Behavior of P/T 4 as a function of 
T/T c . The quasi-particle pressure is obtained by using the 
fitting parameters for z g and z q listed in Table I. Lattice Data 
are also shown as points. 



FIG. 4. (Color online) Behavior of e/T 4 as a function of T/T c 
in the quasi-particle model. The lattice results are shown 
as points and the solid line shows the quasi-particle result 
obtained by utilizing the fitting parameters, a 9>? ,6 9i9 , and 
a!g tq ,b' g q listed in Table. I. 



Moreover, these effective fugacities can be expressed in 
terms of effective quasi-particle number densities. These 
number densities leads to a simple Virial expansion for 
the EOS which we shall discuss in the next section. 

Next, we look at the group velocity of quasi-partons. 
The group velocity can be obtained as v p = dpw p . It 
is easy to see that the modified term in the dispersion 
relations is purely temperature dependent, therefore it 
will not change the group velocity of a quasi-parton. 
v p = p for quasi-gluons and quasi-quarks (u, d) and 

v p = -, p for strange quarks. The dispersion re- 
VP 2 + m2 

lation in Eq.(|S]) contributes to the trace anomaly in hot 
QCD which we shall discuss soon. 

Let us now discuss the significance of of gluon con- 
densate, in the hot QCD thermodynamics. In this con- 
text, D'Elia, Giacomo, and Meggiolaro [f3(| have stud- 
ied the electric and magnetic contributions to the con- 
densate, and shown that near T c the former vanishes, 
however the latter remains unchanged. These authors 
investigated such effects by analyzing the two-point cor- 
relation functions both in pure-gauge sector, and the full 
QCD [ini ■ It has been shown in [52[ that the effects 
of the gluon condensate are significant for T > T c , and 
becomes vanishingly small beyond 2T C . Therefore, in the 
effective mass description of hot QCD for the temper- 
atures, T = 1 — 2 T c , one needs to consider the con- 
tributions of the condensate, while comparing the pre- 
dictions on thermodynamic observables with the lattice 
QCD data. However, in our model, z g , z q capture these 
effects, and we do not need to incorporate the contri- 
butions separately. This can be understood as follows, 
In the lattice data employed here, normalization of the 
pressure and energy density were chosen such that at 
T = 0, these quantities vanish and the effects of 
the gluon condensate may be significant at higher tem- 



peratures. These effects are well captured in the trace 
anomaly in lattice QCD, which is the basic quantity com- 
puted in the lattice. And, all other thermodynamic quan- 
tities have been derived from the trace anomaly. In turn, 
these effects are automatically encoded in the pressure, 
the energy-density etc.. In our study, such effects have 
been captured in the effective fugacities, z g / q from the 
beginning, since, we have determined them from the lat- 
tice data on pressure. The gluon condensate contribute 
significantly to the energy density, entropy density, and 
the trace anomaly through the temperature derivatives 
of the Zg/g, in terms of modified dispersion relations. 



B. Viability of the model 

As it has been already emphasized in the previous 
section that the model yields (2+l)-flavor lattice QCD 
EOS almost perfectly. To check further the viability 
of the model, we study the temperature dependence of 
the quasi-particle pressure, energy density, and the trace 
anomaly and check them against the direct lattice re- 
sults. Lets us first discuss the temperature dependence 
of the pressure. We have plotted the quasi-particle pres- 
sure along with the lattice data in Fig. 3. We find that 
the agreement between the lattice data and quasi-particle 
model for the EOS is almost perfect beyond T c . The tem- 
perature dependence of the energy density as a function 
of temperature is shown in Fig. 4. The quasi-particle 
results agree well with the lattice data beyond T c . 

Encouraged from the crucial observation that lattice 
and quasi-particle model predictions are the same, we 
now proceed to study the trace anomaly as a function of 
temperature obtained from the quasi-particle model. 



6 




FIG. 5. (Color online) Behavior of n g as a function of T/T c . 
N g is obtained by employing the discrete data points for z g 
as well as fitting parameter for of z g listed in Table. I. The 
Stefan-Boltzmann (SB) limit of n g is also shown. 



1. The trace anomaly 



exactly by the numerical evolutions of the integrals in 
Ea. ([T5]) . They are represented in terms PolyLog func- 
tions merely to understand the non-trivial temperature 
dependence of the effective number densities. Here, we 
see that n g /T 3 , and n q / s /T 3 scales with T/T c in a non- 
trivial way. Their behaviors with temperature are shown 
in Figs. \E\ and [6] 

We have also shown the Stefan-Boltzmann value 
for the number densities in Figs. [5] and |H1 
(n g /T 3 )\ SB = ^C(3)/tt 2 ; (n q /T 3 )\ S B = 3^C(3)/4^ 2 ; 
and (v s /T 3 )\ SB = 3j/ s C(3)/4tt 2 . Note that the SB limit 
has not been obtained even at T = 5T C . The effective 
number densities are roughly ~ 10% away at this tem- 
perature in all the three sectors. This is just the reflection 
of the fact that lattice EOS itself is away from SB limit 
there. 

Using Table. I, the quantities Tdr ln(%) and Td\n(z q ) 
can easily be obtained as, 

Td T Hz g ) 
Td T Hz q ) 



5b g 
2b' 



X< X„ 



x> x„ 



~ \ 2b' 



X< X q 
X> X a 



(14) 



Trace anomaly gets contribution from all the three sec- 
tors. The gluonic and light quarks contributions come 
purely from the modified part of the dispersion rela- 
tions. On the other hand, in the strange quark sec- 
tor trace anomaly gets additional contribution from the 
mass. We denote the trace anomaly by A = (e — 3P) = 



Td T \n{z g ) 



T 3 



A ff . , , n„ 

^L=Td T Hz q )^L 

J± =Td T HZg)^ 



(12) 



where n g , n q and n s are the effective number densities 
for the quasi-partons and arc defined by, 



2tt 2 

Vq_ 

2tT 2 



dpp 2 ff 



"I 



T 



is g PolyLog[3, z g 



dpp 2 f q eq = T A 



2 rs 



^2 I terror 



-v q PolyLog[3,-z q ] 

IT 2 

-is s PolyLog[3, ~z q ] 



3m 



> > 2 



ln(l + z q ). 



(13) 



We shall first discuss the behavior of these effective num- 
ber densities as a function of temperature and thereby 
the temperature dependence of the trace anomaly in 
hot QCD. We determine the effective number densities 



It is now straight forward to compute the trace 
anomaly employing the fitting functions for the effective 
fugacities; and effective quasi-particle number densities 
using Ea.(|12p. Behavior of A along with corresponding 
lattice values has been shown in Fig. 7. As it is clear 
from Fig. [3] and Fig. [71 pressure and trace anomaly 
computed by employing the quasi-particle model show 
good agreement with the lattice data of the same. This 
sets the utility of the model. Once these two quanti- 
ties are known, it is straight forward to determine the 
energy density (e) and the entropy density(s) by the 
standard thermodynamic relations, e = 3P quas i + A, 



AP -4- A 

^ J quasi i 



If we see, closely the behavior of energy-density or the 
trace anomaly as a function of temperature in Figs. 4, 
and 7, we observe that around T = 1.7T C , both the quan- 
tities are not showing smooth behavior. There is no phys- 
ical reason associated with this. This is merely the ar- 
tifact of the two distinct fitting functions for z g and z q 
below and above this temperature. At this temperature 
these functions have different slopes. This problem may 
not be present, if we could have found a single fitting 
functions for z g and z q for the whole range of tempera- 
tures. As emphasized earlier, these fitting functions were 
needed to compute, the temperature derivatives of z g and 
z q . Moreover, these points merely reflect the fact that at 
this point, perhaps various thermodynamic quantities are 
changing their slopes, but smoothly, which is not cap- 
tured appropriately in the fitting functions. This may 
not be thought of as a serious problem since, one could 
compute the temperature derivatives of z g and z q , di- 
rectly by inverting relations in Eq. (12). In that case, 
we shall not get a continuous curve for energy density, 
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T/T c T/T c 

FIG. 6. (Color online) Behavior of n q (left) and n s (right) as a function of T/T c . They are obtained by employing the discrete 
data points for z q as well as fitting parameter for of z q listed in Table. I including the strange quark mass dependence. The 
Stefan-Boltzmann (SB) value is also shown beyond 3T C . 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



T/T c 

FIG. 7. (Color online) Behavior of trace anomaly as a function 
of T/T c . The solid line denotes the values obtained from the 
quasi-particle description and the points denotes the lattice 



rather discrete points and they will perfectly match with 
the lattice predictions. 

IV. IMPLICATIONS OF THE MODEL 

A. Virial expansion for hot QCD 

To translate QCD interactions in RHIC era in terms 
of a Virial expansion is quite a non-trivial task. There 
are very few attempts in this direction (54[. Here, we 
see that the quasi-particle understanding of equation of 
state for QGP plays very crucial role to obtain a very 
simple Virial expansion of the EOS in terms of effec- 
tive quasi-particle number densities. The model tells us 



that the Virial expansion in (2+l)-navor QCD could be 
subdivided in three sectors, viz., gluonic, light-quark, and 
strange quark sector and one can define the Virial expan- 
sion in each sector and finally combine them. 

We begin with the purely gluonic sector first and sub- 
sequently discuss the matter sector (light quarks and 
strange quarks). 

1. Virial expansion in purely gluonic sector 

To obtain the Virial expansion in this sector, we Taylor 
expand the pressure, P g and the effective gluon number 
density, n g in gluonic sector in the powers of z g (assum- 
ing z g < 1) and eliminate the explicit dependence of z g 
from both the expressions. This technique which is the 
standard way to obtain the Virial expansion of a ideal 
Bose/Fermi gas with fixed number of particles [Hj], is 
equally applicable here, although the systems are phys- 
ically distinct, since z g and z q do not correspond to the 
particle conservation. The procedure is straight-forward. 
For the sake of completeness, we shall write a few steps. 
The expressions for P g and n g in the powers of z g are 
obtained as, 

,15) 

i = ^X>4. tie) 

9 th j =1 

The quantity Xth = 1/T is the thermal wave length of 
gluons and the coefficient bi is given as, 

1 f°° 1 

bl = 2^1J ««pH«)** s ^4- ( 17 ) 
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Consequently the ratio of Eq. (fT5"j) with Eq. p^]) and 
expanding it as, 



n g T 



k=l 



where iV 9 = 2a and afc's are the Virial coefficients. Using, 
Eq. (fTSjl and Eq. lfTB]) in Eq. (fTS|) and comparing the terms 



of order z g , z 2 and z 3 , we obtain, 



P 1 00 
u,T A 3 ^ v y * 



th /_ 



00 



A 3 



1 OO 

i-£(-l)'-iJ6{4. 



i=l 



(22) 

(23) 
(24) 



62 _ ^ 2 



a 3 = 



26, + 



"4 



6? 

1 ( , 206 3 



186 2 6 3 



(19) 



One can, in principle obtain all the Virial coefficients 



comparing various order coefficients of z g in Eq. ([18l) 
Now, the Virial expansion up to 0((N' A 3 ^) 3 ) can be writ 



ten as, 



T 



1 + a 2 {N'X h ) + a 3 (N' A 3 ,) 2 + a 4 (N' A 3 ,) 3 . (20) 



If we exploit the temperature dependence of effective 
gluon number density shown in Fig. 5. We can compare 
the strength of various order terms in Ea. ((20|) . We find 
that third term is << second term and fourth term is 
<< third term. Remember that the Virial coefficients in 
Eq. (|19[) do not acquire any temperature dependence and 
are same as those for an ideal gluonic plasma with tem- 
perature independent fugacity. The interactions merely 
renormalizes the number density of quasi-gluons. This 
confirms our view point that hot QCD medium effects 
can entirely be mapped in to the non-intcracting/wcakly 
interacting quasi-particlc degrees of freedom. The valid- 
ity of the Virial expansion in this sector is ensured by the 
fact that N' g Xf h « 1. 

Let us now move to the matter sector and first dis- 
cuss the light quark sector followed by the strange quark 
sector. 



2. Virial expansion in the matter sector 

We shall exactly follow the same procedure discussed 
earlier to obtain the Virial expansion. We denote the 
light quarks contribution to pressure as P q and contribu- 
tion from strange quarks to P s . Expanding these quanti- 
ties along with effective number densities n q and n s in the 
power of z q ( assuming z q < 1), we obtain the following 
expressions, 



A. 

VqT 



1-1 



h z 



(21) 



Where the coefficients 6! is defined as, 



6 -2^7 



u 2 exp 



rh 2 )du 



(25) 



Repeating the analysis same as for gluons, and denot- 
ing the Virial coefficients in light quark sector as aj and 
strange quark sector as a ; s , the Virial expansion for P q 
and P s would have the following forms, 

A = a\ + al{N'\\ h ) + 4(N'X 3 th ) 2 + aWA 3 ,) 3 



n q T 

Ps 



= a s 1+ a s 2 (N' s X 3 th ) + a s 3 (N' s \* th y + ai(N' s \° h ) 



1 \3 \2 



jl \3 \3 



(26) 



where N' q 



'f- and N' s = ^f-. In principle, all the Virial 

coefficients (af and af) are possible to compute in terms 
of bi and b\ . We shall only discuss up to the fourth Virial 
coefficient. These coefficients are obtained as follows, 



q T q b 2 7T 5 

i? = 1; ak = -77 = — 



-2&3 + f 



1 "-3 



364 



b\ 2 4 
206 3 , 186 2 6 2 



bf 



b\ h 



a\ = 1; o$ = (b 2 — )/(& 



2 5 7T 2 



2 2 7T 2; ' 



2b' 



(b[) 4 



(b[) 3 

, , 20(6 2 ) 3 18^6!, 



(27) 



Again the dominant contribution is from the second 
terms in the Virial expansion of P q and P s in Eq. (|26|) .This 
we have observed by exploiting the temperature depen- 
dence of the effective quasi-particlc number densities 
shown in Fig. 6. 

Since the total quasi-particle pressure is P = P g + P q + 
P m , Virial expansion of the full EOS can be obtained by 
using the individual Virial expansions obtained in the 
effective gluonic sector (EGS), and the matter sector ( 
Eas.([2"D")l and ([2l)]>. Note that all the Virial coefficients 
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in EGS as well as in the matter sector are independent 
of temperature. The information about the interaction 
has been captured in the effective quasi-parton number- 
densities. The second Virial coefficient is negative in the 
gluonic sector and positive in the quark sector. This 
is expected from the quantum statistics of quasi-gluons 
and quasi-quarks. It is straightforward an exercise to 
determine the other thermodynamic obscrvables in term 
of effective quasi-parton number densities by using the 
well known thermodynamic relations. These expressions 
will also contain the temperature derivatives of effective 
number densities in addition. In other words, both z g 
and the modification factor to the dispersion relations, 
Tdrhi(Zg), T9rln(z 9 ) will appear in their expressions. 
Finally, the validity of the Virial expansion in the matter 
sector is ensured by the fact that N' q s X^ h << I. 

This is perhaps the first time, we have obtained such 
a simple Virial expansion for hot QCD where interac- 
tions appear as suppression factors through the effective 
quasi-parton number densities. This has only been pos- 
sible due to the quasi-particle description of hot QCD. 
Interesting enough, such a description works well down 
to temperatures which are of the order I.0T c . The Virial 
expansion here highlights the role of interactions in hot 
QCD. The Virial expansion may possibly play crucial role 
to explore a quantitative understanding of Fermi liquid 
like picture of hot QCD interactions as indicated by our 
quasi-particle description and also play important role 
to develop effective field theory and effective kinetic the- 
ory for such a quasi-particle model. We shall leave these 
interesting issues for the future investigations. 

At this juncture, we wish to mention that there has 
been a very recent attempt [55| to study the nuclear mat- 
ter EOS at sub-nuclear density in a Virial expansion of a 
non-ideal gas. The Virial expansion is obtained by con- 
sidering the fugacitics for various species such as neutron, 
proton etc. The method to obtain the Virial coefficients is 
standard one as employed in the present work. However, 
the major difference between the two is in the physical 
meaning of the fugacities. 



B. Comparison with other approaches 

We now intend to compare our quasi-particle model 
with other existing models. In the recent past j40l - l43| and 
in a very recent work [44j . effects of hot QCD medium 
have been interpreted in terms of single particle states 
(effective gluons/quark-anti-quarks) via the Polyakov 
loop. In these approaches, the expectation value of 
the Polyakov loop appears in the effective gluon/quark- 
antiquark distribution functions [43| . It provides a sup- 
pression factor in the form of a effective fugacity to an iso- 
lated particle with color quantum numbers. On the other 
hand, there have been successful attempts to encode the 
high temperature QCD medium effects in terms of ef- 
fective thermal masses for quasi-partons [3514391 ]. In the 
recent past, effective mass models have been employed 



to describe (2 + l)-flavor QCD [H, [57| and the agree- 
ment was found to be good. The effective mass models, 
which we discussed so far are based on the lowest order 
results in perturbative QCD, equivalently leading order 
HTL results. These models are improved by incorporat- 
ing the next order HTL contributions by Rebhan and 
Romatschle (58[- Their predictions were shown to be in 
agreement with the lattice results including (2 + I)-flavor 
QCD. Furthermore, there are other approaches which 
also involves quasi-particle picture of hot QCD along with 
the contribution from the gluon condensate [52j. We shall 
compare our model with these approaches one by one. 

Let us consider the Polyakov loop approach first. 
There are certain similarities and a number of differences 
between our model and this approach. The similarities 
are, (i) both the approaches lead to an effective descrip- 
tion of hot QCD in terms of free quasi-particlcs, (ii) the 
expectation value of the Polyakov loop which plays the 
role of effective fugacity as well as effective quasi-parton 
fugacities (z g / q ) in our model appear as the suppres- 
sion factors in the corresponding quasi-parton distribu- 
tion functions, (iii) in both the approaches the group ve- 
locity (v gr — dpiUp), remains unchanged, (iv) both the 
models are quite successful in reproducing the lattice 
data on thermodynamic observables (For more details on 
Polyakov loop method see (42|. Our model yield lattice 
EOS for SU(3) pure gauge theory almost perfectly which 
the deviations which are one part in a million |26[), and 
same is true for the (2+ 1) -flavor lattice EOS, and (v) 
the effective gluon distribution function in [43, EH has 
a similar mathematical structure as our model. In spite 
of these similarities, our model is fundamentally distinct 
from this approach. Our model is purely phenomcnolog- 
ical, and is more in the spirit of Landaus theory of Fermi 
liquids. We list below the major differences between the 
two approaches. 

• The expectation value of the Polyakov loop ap- 
pearing in the single particle distribution function 
does not change the dispersion relation for quasi- 
partons. On the other hand, in our model, we ob- 
tain non-trivial quasi-parton dispersion relations. 
(T 2 d T ln(z g/q ). 

• The Polyakov loop (its phase) appears as an imagi- 
nary chemical potential in the single particle distri- 
bution functions [HI, [44| . This is unlikely to hap- 
pen in our model. The effective fugacities in our 
model cannot be interpreted as chemical potentials 
(real/imaginary) (since there is no conservation of 
particle number). They are introduced merely to 
capture all the interaction effects present in hot 
QCD medium. 

• Employing our quasi-particle model, one can study 
the the bulk and transport properties of hot and 
dense matter (QGP) in RHIC. These studies have 
been reported for pure SU(3) gauge theory in 0, 
|26| and will be presented separately for full lattice 
gauge theory in the near future. 
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Let us now compare our model with the effective mass 
models. In this approach, lattice QCD data for the 
EOS had been interpreted in terms of effective thermal 
gluon mass and effective thermal quark mass. The quasi- 
particle model proposed in this paper is completely dis- 
tinct from this model. The major difference is in the 
philosophy itself. The effective fugacities are not the ef- 
fective masses and they can be interpreted as effective 
mass in some limiting case (p « T 2 8t ln(z g / q )). More- 
over, our approach explores the Fermi liquid like picture 
of hot QCD. Another major difference in two of the ap- 
proaches can be realized in terms of group velocity v gr , 
v gr in two approaches is not the same. In the effective 
mass approaches v gr depends on thermal mass parameter 
(v gr = p/ \/p 2 + m(T) 2 ). We have obtained a Virial ex- 
pansion for hot QCD in terms of the the quasi-partons. 
This has not been done employing either of these two 
models. 

Let us now discuss the quasi-particle models which 
incorporate the effects of gluon condensate explicitly 
in the analysis [lU, [59|. In these studies, the impor- 
tance of the gluon condensate is highlighted, and its 
effect on the thermodynamic obscrvablcs in hot QCD 
was studied in detail [521 ] . In particular, Castorina and 
Mannarelli [52| , have analyzed the thermodynamic prop- 
erties of hot QCD between 1 — 2 T c by explicitly incor- 
porating the gluon condensate along with the gluon, and 
quark quasi-particlcs with thermal masses. The results 
show excellent agreement with the lattice predictions, 
both in the pure glue sector and the full QCD sector. 
Apart from the differences in the dispersion relations, 
and the philosophy with our model, there has been a very 
crucial difference. As emphasized earlier, in our model 
the effect of gluon condensate has been incorporated from 
the beginning, and not treated explicitly as in [52j . How- 
ever, both the models are equally successful to describe 
the lattice QCD thermodynamics. 

There is an alternate way to interpret the effective fu- 
gacities in terms of effective mass, as follows. Let us 
suppose, z g/q = exp(—m e ff\ g / q /T). The quantity, m e // 
can be thought of as m e ff = g'(T)T, where g' is an ef- 



fective coupling. It is to be observed that 



are of 



the order of 0.15 around T c ; it leads to an estimate for 
g' ~ 2.0. This indicates the non-perturbative nature of 
hot QCD matter near T c . Moreover, g' becomes less than 
one beyond 1.3T C , and this observation is valid in both 
effective gluon and matter sector. 



bye mass in terms of the quasi-parton distribution func- 
tions, which are obtained from the (2+l)-latticc QCD 
EOS, is given by, 



Ml 



-2N C Q 2 



d 3 p 



d P f 9 eq + Q 7 



. (28) 

where Q 2 is the effective coupling which appears in the 
transport equation. If one assumes QGP as an ideal sys- 
tem of massless gluons and quarks, Ea. (|28|) reproduces 
the leading order HTL result for the Debye mass; with 
the identification that Q 2 = g 2 (T) (g(T) is QCD running 
coupling constant at finite temperature). 

Employing the distribution functions displayed in 
Eq.([T]) to Eq. (|28p. we obtain the following expressions 
for the 



Mn = Q 2 T 2 



N c 6 Polylog[2, z g 



3 7T 2 

1 -12 PolyLog [2, -z q ] m 2 

2 X V 2 4^2 



ln(l + z q ) 



(29) 

While determining the Debye mass in Eq. ((2T)|) from 
Eq. ([2"5)l , we employ Eq.([T]) and the following standard 
integrals, 

P 2 dp * 9/q ~f P L. 2 = ±^PolyLog[2 1 ±z g/q ] 



z q exp(-/3p) 



Pdp (l + l q cx P R3p))2- 03 



-ln(l + z,). (30) 



The Debye mass with the Ideal EOS(z 9 = 1, z q = 1) 
will be, 

(A/£) 2 =Q 2 T 2 (^ + i-gln(2)). (31) 

To analyze the role of interactions, we define the effec- 
tive charges Q g , Q q , and Q s as, 



Qg=Q 



Qq = Q 



f6PolyLog[2,z g ^ ^ 1/2 



-12PolyLog[2,-z q }^ 1/2 



(32) 



Debye mass could be written in terms of these effective 
charges as, 



Debye screening mass and charge 
renormalization 



To investigate how the partonic charges modify in the 
presence of hot QCD medium, we consider the expression 
for the Debye mass derived in semi-classical transport 
theory (60j in terms of equilibrium parton distribution 
functions. The same expression was obtained from the 
chromo-electric response functions of QGP @]. The De- 



Ml = { f Q 2 + \ (Q 2 + Q 2 ) |r 2 - Q 2 T 2 £j ln(l + z q ). 

(33) 

As stated earlier, in the ideal limit Ea. (j33"|) will reduce 
to Ea. ([5T|) . The quantities, Q g , and Q q ^ s approach to Q. 
This observation tells us that interactions merely renor- 
malize the effective partonic charges. In fact, the effective 
charges are reduced as compared to Q and asymptotically 
approach to the ideal value, Q. 
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FIG. 8. (Color online) Behavior of /j,d as a function of T/T c . 
Note that fid approaches to the ideal value only asymptoti- 
cally. 



To see, how interactions modify Debye screening mass, 
we consider the ratio fid = Mp/M^. The behavior oifid 
as a function of temperature is shown in Fig. [5J fid ap- 
proaches the ideal value unity only asymptotically and 
fid < 1- This implies that the presence of interactions 
suppresses the Debye mass as compared to its ideal coun- 
terpart. 

Next, we compare Mn with the Debye mass obtained 
in lattice QCD. Let us first discuss the Debye screening 
mass computed in lattice gauge theory. It has been cal- 
culated in pure gau ge theory (Nf = 0) [6l|, in 2-flavor 
QCD (Nf = 2) 0, II, and in 2 + 1-flavor QCD [H. 
The lattice data on Debye mass have been fitted with 
the simple ansatz motivated by leading order result on 
Debye mass, to^ = A mjy > , where m 1 ^ denotes the lattice 

data and, m 1 ff = v/(l-i- ~^f)g(T)T denotes the leading 
order Debye mass. Here, g(T) is the two loop running 
coupling constant. This form fits the data quite well if 
Aw 1.4- l^Q. 

In our case (see Eg 155)1 . Q 2 is a free parameter. We can 
fix it to match the Debye mass with the lattice result, 
m^. This leads to, 

C 2 A 2 9 2 (T) 

H {QPolyLog[2,z g ])/^ 

Q 2 = AV(T){((6 PolyLog[2,z g ] 

\ m 2 )~ 1 
-12 PolyLog[2, -z q ])/v 2 \ - — ln(l + z q ) j 

(34) 

in EGS and full QCD with Nf flavors respectively. In 
other words, the Debye mass obtained from the quasi- 
particle model can exactly be matched with the lattice 
results for the Debye mass. The Debye mass is needed 



to determine the transport parameters for quark-gluon 
plasma in RHIC (2(| . It is of interest to derive the form 
of heavy quark potential [&5L [66| employing the formal- 
ism of chromo-electric response functions Q, [65| . While 
deriving the potential, one should keep in mind the fact 
that hadronic phase to quark-gluon plasma transition is 
a crossover [63, |(38[ rather than a true phase transition. 
These issues will be taken up in a separate communica- 
tion in near future. 



V. CONCLUSIONS AND FUTURE PROSPECTS 

In conclusion, a quasi-particle model for (2 + 1)- 
flavor lattice QCD has been proposed which is valid 
in the deconfined phase of QCD. The interactions have 
been encoded in to the effective gluon and quark fu- 
gacities. These effective fugacities non-trivially modify 
the single quasi-parton energies and lead to the trace 
anomaly in hot QCD. The description accurately repro- 
duce the lattice QCD pressure, energy-density, and the 
trace anomaly. In particular, the model accurately repro- 
duce their low and high temperature behavior. We find 
that the model is fundamentally distinct from the other 
quasi-particle models (effective thermal mass, Polyakov 
loop models, and models with gluon condensate). 

Employing the model, temperature dependence of the 
effective quasi-particle number densities has been ob- 
tained. A Virial expansion for QGP has been proposed 
in terms of effective quasi-particle number densities. The 
Virial expansion of the quasi-particle equation of state 
gets contribution from three sectors, viz., the effective 
gluonic sector, the light quark-sector, and the strange 
quark sector. These sectors were dealt separately and 
eventually lead to the complete Virial expansion. This 
is perhaps the first time such an Virial expansion has 
been proposed for hot QCD. Interestingly, the Virial ex- 
pansions came out to be mathematically similar as that 
for an ideal system of gluons, light quarks and, strange 
quark with temperature dependent fugacities. The Virial 
expansion has ensured that the quasi-particle are non- 
interacting. The interactions merely modulate the quasi- 
particle number densities and modify the single quasi- 
particle energies in a non-trivial way. 

The Virial expansion may play important role to ex- 
plore the Fermi liquid like picture of hot QCD in the 
matter sector and, in building effective kinetic theory 
with the quasi-particle model. The Virial expansion has 
revealed that the interactions appear to various observ- 
ables determined by employing the quasi-particle descrip- 
tion only in two ways, either through the effective fugac- 
ities (act as modulation factors) or through the modi- 
fied dispersion relation. Finally, Debye mass has been 
obtained employing the expression obtained from semi- 
classical transport theory and effective coupling has been 
determined in terms of effective fugacities. This observa- 
tion will be required in determining the transport coef- 
ficients for quark-gluon plasma. We find the the Debye 
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mass obtained from the quasi-particlc model can exactly 
be matched with the lattice results. 

The implications of model to study the transport coef- 
ficients (shear viscosity, bulk viscosity) will be taken up 
in the near future. It would also be of interest to ex- 
tend the present model in the case of finite baryon den- 
sity, and studying quark-number susceptibilities. It is to 
be of great interest to establish possible connections of 
our quasi-particle model with the Polyakov loop models, 



which is a matter of future investigations. 
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